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Abstract: We study the sub-Riemannian geodesic problem in the 2n + 1- 
Heisenberg group. An independent set of 2n left-invariant vector fields determines 
a SR structure. We are interested in sub-Riemannian length minimizers for this 
structure. We use the hamiltonian formalism and apply Pontryagin maximum 
principle to write the equations for the extremals and describe some properties 
of the extremal curves. We consider geodesies in the group with a left-invariant 
Riemannian metric. We obtain their equations and prove that the set of directions 
of all rays issuing form zero is the sphere of unit horizontal directions. 
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1. INTRODUCTION lowing optimal control problem: 

n 

Let L be the 2n-M-dimensional nilpotent Lie alge- 7 = X! + ViYi{-f)) (1) 
bra with the following commutation rules in some 

basis Xi,Yi,X2,Y2,...,X„,Yn,T: 7 G G,Ui,Vi G R, 



[X,,Xj] = [Y,Y,] = adT = 0, 

and G be the connected, simply connected Lie 
group with the Lie algebra L. We consider 
Xi,Yi,T as left-invariant vector fields on G. 
The set {Xi,Yi,--- ,Xn,Yn} determines a left- 
invariant sub-Riemannian structure on G: 

D =. span (Xi, Fi, • • ■ , X„, r„), {X,,Y,) 

= {x,,Xj) = (y„ r,) = i, J = 1, 2, • • • , n, 

called the fiat (2n, 2n + 1) sub-Riemannian struc- 
ture. Such structure is unique, up to isomor- 
phism of Lie groups. We are interested in 
sub-Riemannian length minimizers for this sub- 
Riemannian structure, i.e., in solutions to the fol- 



7(0) = go, 7(T) -91 fixed, (2) 




Problem l(I|)-ll3]) was considered by Brock- 
ett(1981), and Liu and Sussmann (1995) for 
n — 1. G is a 2-step nilpotent Lie group and 
L = D (B [L, L]. In this note we consider also 
geodesic lines in the Lie group G with a left- 
invariant Riemannian metric. 

The fiat (2n, 2n + 1) sub-Riemannian structure 
gives a local nilpotent approximation for an arbi- 
trary sub-Riemannian structure with growth vec- 
tor (2n,2n-|- 1), see Montgomery (2002). The dy- 
namics of a classical electric charge in the plan un- 
der the infiuence of a perpendicular magnetic field 



can be described by means of flat sub-Riemannian 
structure. 

In general, if G is a 2-step nilpotent Lie group and 
{D, (,)) is a left-invariant sub-Riemannian struc- 
ture such that L = D (B [L, L], any minimizer is 
normal in some sub-group of G and hence any min- 
imizer is smooth, see Abib (2002). This extends 
known results on the so-called Gaveau-Brockett 
problem, Brockett (1981) and Gaveau (1977). 

2. MODEL 

We choose the following model for the flat 
(2n, 2n + 1) sub-Riemannian structure: 

G — R^"+^ with coordinates q — (x, z) where 

X = (a;i,yi, ,a;„,j/„),z e R, and group law 

deflned as follows: 

{x, z){x' , z') — (x + x' ,z + z' + 

n n 
i=l 1=1 

d d 

= "5 — 

axi az 

oyi dz 
D = KeTe 

n 

6 = dz + ^^{ytdxi - Xidyi) 

i=l 

Control system |[T| has the full rank and the 
state space G is connected, thus the systems is 
globally controllable on G. A standard existence 
theorem from optimal control theory, see Liu and 
Sussmann (1995), implies that any point can be 
connected with qo € G hy an admissible length 
minimizer. 



3. EXTREMALS 

A admissible curve for D is an absolutely contin- 
uous curve on G which is tangent to D almost 
everywhere; applying the Pontryagin maximum 
principle each minimizer parametrized by arc- 
length is a normal extremal or abnormal extremal, 
see Liu and Sussmann (1995). 

3.1 Abnormal extremals 

An abnormal extremal is a admissible curve for 
D which is the projection onto G of characteristic 
curve (in the symplectic sense) of the annihilator 
D-^ of D in T*G. Introduce the linear Hamiltoni- 
ans corresponding to the basis flelds: 

/i,(A) = x(x,), h{x) - A(ro, h{x) = A(r), 



A e T*G. 

is deflned by equations hi = ki = 0, i = 
1, • • • ,n. Assume T{t) G T*G — {0} caracteristic 
curve of D^, then t{t) G span {hi{T{t)), ki{r{t) ; 
i = 1, 2, • • • , n} and h^{r{t)) = h{r{t)) = 0. 
Abnormal extremals are exactly the constant 
curves. 

3.2 Normal extremals 

Normal biextremals are trajectories of the Hamil- 
tonian system 

A = i?(A), A e T*G (4) 

with the sub-Riemannian Hamiltonian H = 
+ ^i)- 111 the coordinates (g, A) ~ 

[xi,''' jXmVl,''' i?/n;^;Cli''' 7Cn;'7lj''' iVnjC) 

of T*G system ^ reads 



ii = £.i + VtC 
yi=T]i- XiC 



(5) 
(6) 

(7) 



z = ^(ii2/i - yiXi) 

i=l 

6 = mC (8) 

Vi^XiC (9) 

C = (10) 

If the constant C is 0, the curve {x{t),y{t)) is 
the straight line, along which the motion takes 
place with velocity constant. The function z{t) 
is then computed by solving ([7]) and the curve 
{x{t), y{t),z{t)) is a straight Hue in G. If C is non 
zero constant, the equations |[8l)-|[9l) imply 



Xi - 2yiC, = 0, 



Hi + 2xiC = 0. 



So 



Xi = Ai cos 2C,t -\- Bi sin 2C,t, 
fji = —Ai sin 2(t + Bi cos 2Ct, 
Ai, Bi 

are constant. The condition that our curve q(t) = 
{x{t), y{t),z{t)) is parametrized by arc-length then 
says that + Bf) = 1. We obtain 

x^{t) - ^(cos6i, -cos(2Ct + 6i,)), 

Vtit) ^ ^ism{2Ct + 9,)-sin9,) 

with Ai = ri sin 9i,Bi = r i cos 9i. The equation 
im imply 



and then z{t) = ^rf{— - 



21 t sin2Ci 



2C (2C)2 



The Pontryagin maximum principle, applied to 
this case, impHes that all the minimizers curve 
parametrized by arc-length are normal extremal 
q{t). 

4. CONTACT TRANSFORMATIONS 

For n = 1, the Lie algebra of symmetries of the 
flat (2,3) distribution D is parametrized by ar- 
bitrary smooth functions of three variables and 
the Lie algebra of symmetries of the flat (2,3) 
sub-Riemannian structure is 4-dimensional, see 
Sachkov (1998). 

For n > 1, let C be the sheaf of germs of inflnites- 
imal automorphisms X such that Lx6 = f6, the 
function / depending on X. C is the sheaf of 
contact transformations. 

Let X ^ C vanishing at 0, and 4>t be the one 
parameter group generated by X. Then 0t(O) = 
and (j)f9 = ftO. Let Vi be the subspace of L 
deflned by 6*0 =0. Then 0t leaves Vi invariant. 
Furthermore (j)^d9 = dft A 9 + ftd9 . This implies 
that the form d9 restricted to Vi is preserved up 
to scalar factor. If we pass to the inflnitesimal 
we see that the linear isotropy algebra of C is a 
subalgebra of the algebra g deflned as follows: 

let w e L not lying in Vi and g = K + N + M 
where K is the set of all A e End(L) with Au = 
and A/Vi £ sp(Fi), 

N is the of all B of the form Bui = for vi e Vi 
and Bu e Fi, 

M is the set of all multiples of the linear trans- 
formation C where Cvi = vi for vi € V\ and 
Cu = 2u. 

The space g^^^ C Hom {g, L) called the flrst pro- 
longation of g, is the set of all T e Hom {g, L) 
which satisfy (T(w))(u) = {T{v)){uj) for all w, 
V d L. In the case, g*-^^ = sp (T^i)*-^-* -t- g and 



sp {Vi 



g'^^ ^) for all k. The algebra 



sp (Vi) is of inflnite type because sp {ViY^^ can be 
identifled with S''^+^(V\*); then g is of inflnite type. 

Corresponding to the coordinates z, xi , • • ■ ,Xn,yi, 
■■■ ,yn, let us choose the basis u — 
eo, ei, • • • , e„, e„+i, • • • , 62 of L. Then it is easy 
to verify that the following vector flelds are 
all inflnitesimal contact transformations (where 
9 = dz + J^'i^iiVid-^i - ^idyi))- 

d d d d d 

oxi oz oyi oz oz 



E 



^ _d_ 



' dxi 



AiJ-{-nXi 



9yj 



where ^ Ars ® e* G sp{Vi), 

r,s—l 

dn r\ i~\ 

\ ^ a \ - a 



dx^ 



d 



^ d d \ d 



J dyj 



dz 



d 

h Xi 

oyi 



d d \ d 



dyj 



dz 



7) z 



>From a) we see that C is transitive. From 
P) we see that the linear isotropy algebra of C is 
indeed g. A direct computation shows that the 
vector flelds a), 13) and 7) form a Lie algebra. C is 
an non-flat, infinite, transitive LAS (Lie algebra 
sheaf) . 



5. RIEMANNIAN CASE 

We deflne the left-invariant metric on G by 
taking Xi, Yi, T as the orthonormal frame. 
Let V the Riemannian connection of left- 
invariant metric an {eh,h = 1, 2, • • ■ , 2n -|- 
1} = ,X„,r„,T}. Then Ve.e, = 

+ 4 - cjj^)efe where [e^e^] = 

For G the matrice (Ve-)ij" is: 



/ 

-T 



T 







T 

-T 

-^n 



Xi 



Xn 





(11) 



Let C(t) = Er=l (^^^(f)X^{t) + V^{t)Y^{t)) + ^{t)T 

c(t) 

geodesies issuing from G G. Then V^^^j = and 
the table give 

Ui + 2jVi = 0, 
ill - 2^Ui = 0, 
7 = 0. 

Therefore, because the parameter t is natural we 
have X]"=i("i + ^i) + 7^ = 1 and we could take 
7(i) — 7 where the constant 7 € [—1,1] is the 
cosine of the angle between c(0) and the T-axe. 

For 7 7^ 0, Ui{t) — rtCos{2jt + 9i) and Vi{t) = 
Vi cos(27i 4- 9i) where r| = uf +vf. 



In coordinates yi, • ■ ■ , x^, yn, ^) the equations 
for geodesies c{t) — • ■ ■ , yn{t), z{t)) are: 

Xi — Hi — 

n n 

z = 7 + Yl ^« Wa;. W - Wy» W (12) 

i=l i=l 

We have 

ri 

Xi[t) — —— cos 2jt + Qi 
27 

V ■ 

.X 37^-1 6, 

^ ' 27 27 ' 

+—ri sin 27i + c 
27 

for some numbers a^, bi, c which could be defined 
from the initial condition c(0) = 0. 
If 7 = 0, then they are horizontal and satisfy to 
the following: 

x,it) = u,{0)t, y,{t) = v,{0)t, z{t) = 

The non-horizontal geodesies (i.e., with 7 / 0) 
are not rays. Indeed. Let I7I 7^ 1. For every 
geodesic c{t) their Xi, yi coordinates are periodic 
functions. Their projection on z = hyperspace 
is a circle. Take two points on c{t) with difference 
of their parameters equals ^. Then, besides the 
geodesies c{t) their is also the vertical geodesic 
z{s) = (0, • • • , 0, s) connecting them. Its length is 
equal to the difference of z coordinates of consid- 
ering points. But, as easy to check this is always 
strictly less than ^, i.e., the length of the interval 
of the c{t) between them, because t is natural 
one by definition. Therefore, every geodesic c{t) 
with I7I 7^ 1 is not minimal. Vertical geodesic 
z{s) (with I7I — 1) also is not minimal, because 
that the vector field Z{s) = Xicoss + Yi sin is 
parallel along it and the sectional curvatures of G 
in 2-dimensional direction generated by Z{s) and 
i(s) equals 1. This means that the index of every 
interval of z(s) with a length greater than tt is 
positive, and z is not minimal. 

Now let us check that all horizontal geodesies 
(with 7 = 0) are rays. Suppose that the hori- 
zontal geodesic a{t) intersects in some point with 



another geodesic c(s) issuing the same point 0. 
Because horizontal geodesies do not intersect, c(s) 
is non-horizontal one. For c{t) from lfT2|) we see 
that the length Z(c) of the projection curve c(t) — 
(ci(t), • • • , C2n(t), 0), measured in euclidean coor- 
dinates, is strictly less than t, i.e., the length of the 
curve c in the left-invariant metric. Than in the 
submanifold z — considered as euclidean space 
the length of every chord, connecting some points 
of c is less than 1(c). If the geodesic c issuing from 
intersects some horizontal geodesic a{t) also issu- 
ing from 0, then a is the chord of the corresponding 
c. For horizontal geodesies their lengths in the left- 
invariant metric coincide with the usual euclidean 
length in z = 0. The argument above means that 
the length of a is strictly less then that of c. There- 
fore, we have the following statement: the set of 
directions of all rays issuing from in the group 
G is the sphere of unit horizontal directions in the 
point i.e., {w = (wi, • ■ • , W2„,0); = 1}. 
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